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We study the metal-to-insulator transition of the Hubbard model at zero temper- 
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is lower in energy for any interaction in the coexistence region and that the transition 
is of second order. 
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Introduction The correlation induced metal-insulator transition (Mott-Hubbard transi- 
tion) is one of the prime examples in which strong correlations dominate the low-energy 
behavior of a physical system. A theoretical treatment of the problem requires an approach 
which is non-perturbative in the interaction. Recently, new insights into the problem were 
gained using the limit of infinite dimensionality It allows for a mapping of a variety 

of lattice models onto impurity problems in a self-consistent ly determined bath f||,[f] and is 
therefore a natural way to formulate a mean-field theory of itinerant systems. While be- 
ing simpler than the original problem, the resulting mean-field theory remains a formidable 
many-body problem which has to be solved using numerical methods. Recently the Hub- 
bard model has been investigated by several groups using Quantum Monte Carlo (QMC) 
simulations and self-consistent perturbation theory (PT). |5|-f§ While a combination of both 
methods established the existence of a Mott-Hubbard transition at a finite value of the 
interaction U in the paramagnetic phase of the Hubbard model at half-filling, important 
questions regarding the nature of the transition remain unsolved. 

In a previous work |]], the coexistence of metallic and insulating solutions over a finite 
range of values of U has been demonstrated. While the metallic solution disappears con- 
tinuously at a value U c2 , the insulating solution disappears abruptly at a value U c \ < U c2 . 
At finite temperature, the difference between the free energy of the solutions is dominated 
by the entropy term. The large entropy, which is a result of the degeneracy of the ground 
state in the insulating case, made it possible to unambiguously determine the existence of 
a first order transition line close to U c i(T). As the temperature is reduced, the free energy 
approaches the energy, therefore an accurate evaluation of the energy is necessary. Depend- 
ing on which solution is lower in energy two very different scenarios may take place: If 
E Ins < E Met , the transition will be close to U cl and the sudden destruction of the metal- 
lic state implies a first-order transition even at T = 0. On the other hand, in the case 
EMet < Ei ns , the metallic solution continuously merges with the insulating one at U c2 , and 
the quasiparticles display a diverging renormalized mass. || 

While the limit T = cannot be attained by QMC simulations, within the second-order 
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perturbative approach the energies of the two solutions are almost degenerate, making the 
consideration of higher-order corrections necessary. An alternative numerical approach to 
the problem was introduced recently: While the large d mean field equations are functional 
equations for the Green function G(iu n ), an approximation can be obtain by modeling G(iuj n ) 
using a finite number N of parameters, which reduces the functional equations to non-linear 
algebraic equations in N unknowns. Following this idea, two different parameterizations 
were introduced. fllOlJlll Both take advantage of a mapping of the lattice problem onto an 



Anderson impurity model with a self-consistently determined bath. The N parameters that 
model G(iu n ) define the hopping amplitudes and energies of the effective electron orbitals 
of the bath, as will be discussed in detail in next section. The resulting problem can then be 
solved at T = by exact diagonalization of the effective Hamiltonian. This is followed by the 
new determination of the set of parameters, and the procedure is iterated until convergence 
is attained. The method is thus non-perturbative in nature and overcomes the problems of 
both QMC and PT, allowing for an accurate evaluation of the energies at T = 0. 

In this paper we apply this approach to the study of the Hubbard model. We establish 
the coexistence of metallic and insulating solutions over a finite range of the interaction 
parameter U and show that at T = the metallic solution has lower energy than the 
insulating one, implying that the metal-insulator transition in the Hubbard model with 
semicircular density of states is of second order. This justifies a posteriori the relevance 
of the earlier studies || of this quantum critical point which captures the essence of the 
Brinkman-Rice transition. 

Methodology In the limit of infinite dimensionality the Hubbard model, described by the 
Hamiltonian 
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H = ~ (*y + W«rCj> + U ]T(n lT - -)(n a - -), (1) 

<i,j> i 

can be reduced to an effective impurity problem, supplemented by a self-consistency con- 
dition. || As in the previous work we focus on a Bethe lattice of infinite connectivity m, 
which in the non-interacting limit corresponds to a semicircular density of states of width 



At, where the hopping parameter t is rescaled in the usual way as t — > Integrating out 
the degrees of freedom other than the origin, one obtains an effective local action of the form 



SeffM=Ej drdr'cl{r)G Q \r-r')c a {r') 



+U jV(n iT (r)-i)(n a (r)-i). (2) 

In the following we focus on the paramagnetic solution at half filling (// = 0). The self- 
consistency condition then reads Gq 1 {wj v ) = iuj n — t 2 G{iuj n ) where G{iuj n ) = — ^ e tulnT < 
T t c(t)c' (0) >s e ff is the local Green function of the Hubbard model once self-consistency is 
attained. As shown in Ref. f| an action of the same form can be obtained from an Anderson 
impurity model by integrating out the conduction electrons. Note that the self-consistency 
condition implies that the role of the hybridization function is played by the local Green 
function itself. The iterative solution now proceeds as follows: G(iu n ) is modeled by a finite 
set of parameters. In terms of the impurity problem, this represents an effective bath for 
the impurity with a finite number of poles. This effective impurity model is then solved 
by exact diagonalization and a new G{iuj n ) is calculated. A new set of parameters is then 
obtained from G{iuj n ) by approximating it by a function with a number of poles equal to 
the number of sites in the bath (this number is in general smaller than the number of poles 
of G(iuj n )). Note that this represents a further approximation of the method (beyond the 
effective Hamiltonian being finite). The whole process is iterated until convergence of the 
parameters is achieved. 

Exploiting these features, two new similar algorithms were proposed recently |T0| , pT| , 
differing basically in the way the new set of parameters is obtained, that is, how the G(iu n ) 
is parametrized by a smaller number of poles. We will consider both schemes and comment 
on their respective advantages and limitations. 

As mentioned, the number of poles of G(iuj n ) is in general larger that the number of sites 
in the bath, therefore this approximation is an essential ingredient of the scheme. Caffarel 
and Krauth [ID| proposed to obtain the new parameters by a x 2 fit of G(iu n ). Starting with 



an Anderson Hamiltonian of the form 



+U(n c i-~)(n cl -^) (3) 
the self-consistency condition becomes t 2 G(ico) = Ea=i i^hr- We thus have to 

minimize 

X 2 = El^n)-E-^! 2 (4) 

iu) n a=l lL ° n €a 

where we sum over frequencies u n = (2n+l)7rT with small fictitious temperature ( T = .001) 
and large cutoff N^Au ~ 2U, to obtain the new set of parameters V a and e a . Note that this 
Hamiltonian effectively describes an impurity surrounded by a "star" of bath electrons. 

An alternative route was introduced in the context of an extended Hubbard model 
PH . This procedure takes advantage of the fact that the Green function G(z) can be 



decomposed into "particle" and "hole" contributions as G(z) = G > (z) + G < (z) with 

G> ( Z ) =< Q K-(g-£o) ct | > alld G< W =< °\ C K +{ H-E 0) C \° > ■ 

The respective contributions can be obtained from a continued fraction expansion as 

< fo f< \ * .J f >= . <f ° /<lf \Z (5) 
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where |/> >= f\gs >, |/< >= f\gs > and \f n+1 >= H\f n > -a n \f n > -%\f n -i >, 
a n =< f n \H\f n >, b 2 n = </ < ' / '"|/" > 1> ) ^o = 0. This implies that G > and G < can be regarded as 
resulting from a Hamiltonian describing an impurity coupled to two chains with site energies 
a^J < and hopping amplitudes . Again the number of poles in the Green function is in 
general larger than the number of sites of the Hamiltonian and in order to close the self- 
consistency, the continued fraction expansion has to be truncated. The approximation in 
this scheme relies on the fact that the continued fraction representation captures exactly 
the moments of the energy of the Hamiltonian, up to the order retained in the continued 
fraction. It can thus be thought of as a moment by moment fitting. This scheme has 
the numerical advantage that it avoids the multidimensional fit of the Green function, but 
the disadvantage that it can be implemented practically only in the case of a semi-circular 



density of states. In the metallic case an explicit extra site at the Fermi energy is introduced 
in order to better represent the low frequency region and, more importantly, to allow us to 
feed-back a metallic bath. The hopping parameter to this extra site is calculated by a single 
parameter minimization of the expression 



X 2 (a)= J2 \G A (iuj n ,a)-G(tuj)\ 2 (6) 



where now Gvi(iu; n ,a) = ^- + (1 — a)Gisi c {iuj n ). Gn c is the truncated Green function to 
length Nc = Nsite/% and ujl and ujh are low and high energy cut-offs defined by the lowest 
poles of G and Gn c , respectively. In this case the moments will be modified by a small 
factor (a) which decreases as the system size is increased. 
The effective Anderson model therefore reads 

Nc-l 

# = E E ( E <ctcL + b p ,{c p JU + h.c.) 

cr p=>,< a=l 
N c -2 , , 

+ E (K^A+ia + h.c.)) + U(n n - -){n n - -) 

a=l " 

+ Y.Wl^ + h.c). (7) 

In both schemes, groundstate wavefunction and groundstate energy of the Anderson 
Hamiltonian are determined by exact diagonalization (up to six sites) and the modified 
Lanczos technique fl2|]. Systems of up to ten sites can be handled on a workstation. The 



zero temperature Green function of the local site is finally obtained from a continued fraction 
expansion using the recursion method discussed above. 

As mentioned in the introduction, a further advantage of the formulation of the problem 
in terms of an Anderson impurity model is the fact that the energy of the Hubbard model can 
be obtained directly without frequency summations using Anderson model relations. The 
kinetic energy per site of the Hubbard model is given as T = J2<j,k> tGjk{iu) n )e' luJn0+ . 
Taking the limit of infinite coordination number this reduces to T = ^ Yhu„ G(iuj n ) 2 e luln0+ . 
Using the self-consistency condition as well as the the fact that in the Anderson model 

= T,aa v a < flc aa + h.c. > we obtain 



T = y £V«Ile<0\fic aa \Q>, 



(8) 



where a labels the sites neighboring the impurity. The potential energy of the Hubbard 
model is obtained as 



Results In our analysis we have focused on two major aspects: the determination of 
a region where two solutions are allowed, and the resolution of controversy regarding the 
lowest energy solution. The two approaches considered yield a consistent picture of the 
transition. We are able to obtain converged metallic and insulating solutions for a finite 
range of the interaction U within both schemes. We further demonstrate that the metallic 
solution is lower in energy in the whole coexistence region. The energy difference between 
the solutions goes to zero as U C 2 is approached, implying that the transition can be classified 
as second order. This should be contrasted with the results from second-order perturbation 
theory, where the two solutions were found to cross in energy at an intermediate value of 
the interaction U. A point worth noticing (as was already observed within the perturbative 
approach) is that the energy difference between the solutions is much smaller than any 
energy scale of the problem. This is due to an almost perfect compensation of the gain in 
derealization (kinetic) energy, by the loss of energy through double occupancy (potential 
energy), in the metallic state compared to the insulator. 

Metallic and insulating solutions for U = 2.7 inside the coexistence region are respectively 
shown in Fig.|l| and Fig.|2| (the half-bandwidth 2t is set equal to unity). In the first figure 
the Green function displays a narrow resonance at low frequency (note that the pinning 
condition at u = is fulfilled fL3|), while the insulator in the second case merely consists 
of high energy features (upper and lower Hubbard bands). The figures also illustrate the 
consistency of the two schemes considered here. In both, the metallic and insulating, cases 
the agreement is very good. We also find that the results of both methods for the single 
particle Green function on the imaginary axis compare very well with the second-order 
perturbative calculation [|J and QMC |||7| (the latter is discussed in Ref. [|10| ). 



V = U<0\n f ^n fl \0>. 



(9) 
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The kinetic, potential and total energies for the two solutions in the coexistence region 
are displayed in Fig|| An interesting feature is the already mentioned almost perfect can- 
cellation of derealization and double occupancy energy. Another important observation is 
that while a finite size effect is apparent in the results for the kinetic and potential energy, 



the convergence of the total energy is much faster [H|. A few runs for a ten site system 
show almost no difference to the results for eight sites. 

The energy difference of the two solutions is shown in the inset of Fig|J. As the critical 
point U C 2 is approached from below, the finite size effects become relevant for U ~ 2.8. 
This limitation of the scheme is due to the fact that as the low energy scale associated 
with the quasiparticle peak goes to zero close to the transition, the discrete nature of the 
approximation starts playing an important role and the Kondo resonance is represented by 
only a single pole. 

The smallness of the difference in energy between the metal and the insulator can be 
understood from the picture of a second-order critical point where the metallic and insulating 
solutions merge with a vanishing scale A ~ U C 2 — U. The problem can be formulated from 
a variational point of view, with the free energy F becoming an extremum at the metallic 
and insulating solutions, i.e., = jg- = 0. Since the two solutions merge at the point 
Uc2, F can be expanded in power series of Gm — Gj as 

F M -F I = ^(G M -G I ) 2 . (10) 

As the difference between the metallic and the insulating solution is parameterized by A, 
and the second derivative vanishes at the critical point as A, it follows that the energy 
difference goes to zero as A 3 . The critical region cannot be accessed by the present method. 
In order to capture the vanishing energy scale, a higher resolution (i.e. an effective bath 
with more sites) is needed. 

Finally, we would like to comment on the disappearance of the insulating solution at 
U c \. From the "two chain" scheme, the insulating solution is found to persist all the way 
down until the gap closes. This differs from the results of perturbation theory and resembles 
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the Hubbard III scenario for the destruction of the insulating state. [|17],|ll| m the case °f 



the "star configuration", while a converged insulating solution can be obtained at values of 
the interaction U well below U C 2, the question of the closing of the gap cannot be answered 
conclusively. 

Conclusion We have resolved the standing questions regarding the metal-to-insulator 
transition in the Hubbard model in infinite dimensions, using a powerful algorithm to obtain 
Green functions at zero temperature [ 15| . We were able to demonstrate the existence of a 



region in which metallic and insulating solutions coexist, which is in agreement with previous 
results, and showed that the metallic solution is always lower in energy. This implies that 
while at finite temperature the transition is first order, it becomes second-order at T = 0, 
similar to the work of Brinkman and Rice in the context of the Gutzwiller approximation [16| 



||. Since the method presented is very general as well as simple, especially when compared 
to Monte Carlo simulations, it is an appealing approach to the study of strongly correlated 
electron systems. 
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FIGURES 

FIG. 1. Comparison of the metallic and insulating Matsubara Green functions for U = 2.7, as 
obtained from the two variations of the algorithm. Full line "star geometry" and dotted line "two 
chain geometry" (10 sites for the metallic case and 8 sites for the insulating) 

FIG. 2. Kinetic, potential and total energy for the metallic and insulating solutions in the 
coexistence region. Difference between the metallic and insulating solution (inset). From the "two 
chain geometry" . 



12 



